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The search for a new source of CP violation is one of the most important endeavors in particle 
physics. A particularly interesting way to perform this search is to probe the CP phase in the hrr 
coupling, as the phase is currently completely unconstrained by all existing data. Recently, a novel 
variable 0 was proposed for measuring the CP phase in the hrr coupling through the r ± —» ■k 3si 'k 0 v 
decay mode. We examine two crucial questions that the real LHC detectors must face, namely, the 
issue of neutrino reconstruction and the effects of finite detector resolution. For the former, we find 
strong evidence that the collinear approximation is the best for the 0 variable. For the latter, we 
find that the angular resolution is actually not an issue even though the reconstruction of 0 requires 
resolving the highly collimated 7r ±: s and 7r°’s from the r decays. Instead, we find that it is the 
missing transverse energy resolution that significantly limits the LHC reach for measuring the CP 
phase via 0. With the current missing energy resolution, we find that with ~ 1000 fb _1 the CP 
phase hypotheses A = 0° (the standard model value) and A = 90° can be distinguished, at most, 
at the 95% confidence level. 


I. INTRODUCTION 

The Standard Model (SM) of particle physics has been 
extremely successful in explaining all microscopic phe¬ 
nomena observed down to distances of order ~ ITeV -1 . 
However, the model has two notable under-explored cor¬ 
ners: the neutrino and Higgs sectors. Having discovered 
the Higgs boson the next major goal of the Large Hadron 
Collider (LHC) is to scrutinize the properties of the Higgs 
boson and search for deviations from the SM predictions. 

A particularly interesting question is whether CP is a 
good symmetry of the Higgs sector. Since all Higgs boson 
couplings are predicted to be CP even in the SM, any 
observation of a CP-odd or CP-violating component will 
constitute unambiguous evidence of physics beyond the 
SM. CP-violating Higgs interactions may be motivated 
by the fact that baryogenesis requires the existence of a 
new source of CP violation and the possibility that such 
a primordial source of new CP violation may leave traces 
in the couplings of the Higgs boson. 

Recently, a novel variable 0 for measuring the CP 
phase in the h-r-r coupling was proposed in Ref. if]. 
The distribution of 0 was shown to have the simple form 
c — a cos(0 — 2A), where A is the phase of the h-r-r 
Yukawa coupling, with A = 0° and 90° corresponding 
to the CP-even and -odd limits, respectively, while a 
and c are independent of A and 0. What distinguishes 
the h-r-r Yukawa coupling from other couplings is that 
its phase (i.e. the CP angle A) is currently completely 
unconstrained by data [21, even though its magnitude 
(i.e. the rate of r + -r~ production) is roughly consistent 
with the SM prediction USE]. Field theoretically, as dis¬ 
cussed in Ref. [T] , an 0(1) A in the h-r-r coupling can 
easily coexist with other Higgs boson couplings whose 
CP phases are known to be small, such as h-t-t and h-b-b 
(A t < 0.01 and A{, < 0.1, respectively [2]), and h-Z-Z 
(for which Az = 90° is disfavored [SHZ])- In Ref- D3 using 
parton-level analyses it is estimated that the h-r-r CP 


phase, A, may be measured from the 0 distribution to an 
accuracy of ~ 11° at the 14-TeV LHC with ~ 3000 fb^ 1 
of data, or the A = 0° and A = 90° cases may be dis¬ 
tinguished at 5a level with ~ 1000 fb _1 . However, these 
conclusions are for an ideal detector. In this paper, we 
examine the impact of realistic detector effects on the 
measurement of A via 0. 

There are various ways that detector effects can affect 
those estimates. Most importantly, the reconstruction of 
the 0 variable requires resolving the charged and neu¬ 
tral pions in the d —> decay mode, but 

these d are highly boosted as they come from higgs de¬ 
cays. One therefore expects the two photons from each 
7T° decay to hit the same tower of the electromagnetic 
calorimeter (ECAL), thereby making the standard iden¬ 
tification of a 7T° using photons useless. Moreover, one 
expects the same ECAL tower to be hit by the , raising 
the question of how much of the ECAL tower’s activity 
should be attributed to the 7r°. 

Another obvious issue is the reconstruction of the neu¬ 
trinos. In Ref. [J, the so-called collinear approxima¬ 
tion [Sj is adopted, where the neutrino from each r decay 
is taken to be exactly collinear with its parent r. As 
shown in Ref. [I], the collinear approximation results in 
a 60-70% reduction of the amplitude of the cos(0 — 2A) 
curve. So, one might expect that here lies a significant 
opportunity for improvement. In this paper, however, 
we present a theoretical calculation that strongly sug¬ 
gests that it is extremely difficult, if not impossible, to 
outperform the collinear approximation for the purpose 
of 0 reconstruction. 

Adopting the collinear approximation, we then focus 
on the aforementioned issues of pion identifications in 
the presence of realistic detector effects, and propose a 
realistic algorithm to identify “r candidates” specifically 
suited for constructing 0. Surprisingly, we find that the 
amplitude of the cos(0 — 2A) curve of the signal is vir¬ 
tually unchanged by the finite angular resolution of the 
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LHC detectors. Instead, we find that the measurement 
of A through the 0 distribution is significantly ham¬ 
pered by the contamination by the Z —► rr background. 
This degradation arises due to resolution effects in recon¬ 
structing the rf invariant mass as a consequence of large 
uncertainties in the measurement of the missing trans¬ 
verse energy. Fortunately, this is an area where future 
improvement can be hoped for, unlike the angular reso¬ 
lution. We leave the important question of pileup effects 
for future work. 

While we focus on the measurement of A through the 
0 variable, similar but different acoplanarity angles have 
been studied both for e + e^ colliders [SHU] as well as for 
the LHC [T51 - BD] . For the —> 7r ± 7r 0 ^ decay mode, 

which is the most dominant decay mode of the r, the 
acoplanarity angles used in those studies are not as ef¬ 
fective as the 0 variable at e + e _ colliders in which the 
neutrino momenta can be reconstructed up to a two-fold 
ambiguity JT] . An advantage of the LHC studies men¬ 
tioned above is that they can be applied to all 1-prong r 
decays including leptonic decays. However, realistic de¬ 
tector effects are not considered in those studies. Ref. m 
performs a realistic study for the LHC but assuming uni¬ 
versal CP violation for all SM fermions, unlike our study 
in which we assume A is non-zero only for the r. Fi¬ 
nally, the use of the r impact parameter associated with 
a displaced r vertex was studied in Ref. [ 22 ] in the e + e“ 
context. Our 0 variable does not depend on the mea¬ 
surement of the r impact parameter. 

The paper is organized as follows. In Section [TT] we 
review the 0 variable and rewrite it in a simple form by 
expanding it in terms of a small parameter A that char¬ 
acterizes the degree of collinearity of the decay products 
of the taus, with A —> 0 in the exactly collinear limit. We 
use this expansion to argue that the collinear approxi¬ 
mation appears the best we can do for the purpose of 
reconstructing 0. Section m describes the parton level 
simulation of the signal and background processes. Sec¬ 
tion Ev] describes the detector simulation and selection 
of the parametrized detector quantities corresponding to 
the Higgs boson decay products. Section |V| evaluates the 
results of the simulation and sensitivity as a function of 
integrated luminosity, and we conclude in Section |VT| 


II. THEORY 


A. The CP phase A and the variable 0 


Following Ref. pQ, we define the CP phase A by 


C 


r-Yukawa 


Vt 

V2 


hr (cos A + iy 5 sin A)r, 


( 1 ) 


where the scalar field h and Dirac spinor r describe the 
Higgs boson and r lepton, respectively. In order to main¬ 
tain the SM decay rate for /i -> rr, we fix the h-T-r 
Yukawa coupling, y T , to its SM value yf M = m T /v with 


v ~ 174 GeV, but keep the CP phase A as a free param¬ 
eter to be measured. The SM prediction is A = 0. 1 
We are interested in the decay chain 

h -0 T + T~ 

—¥ p + P p~v (2) 

—► 7T + 7 T°P 7T~Tr°iy . 


As calculated in Ref. 1|, the squared matrix element for 
this decay has the form 

| Ad h —tzttv | OC C Q cos(0 2 A), (3) 


where a and c are both positive and depend neither on 
0 nor A. The angular variable 0 is defined through 

cos © a (ki-p T2 ) (k 2 -Prl) - (Pt 1 -Pt 2 ) (ki-k 2 ), 
sin 0 a e^pa fcf p v rX p° 2 , 

with a common proportionality factor. The subscripts 
1 and 2 refer to the r + and t~ branches of the decay, 
respectively (e.g., p^ x = p^ + and p^ 2 = p^_). The 4- 
momenta fcf (i = 1, 2) are defined by 

ki = UiQi + r Pvn ( 5 ) 

with 


and 


u _ a u 

Qi = Pl±i - TC o 


2 Qi-Pri - 4ml 

- T = ——- 

+ m 2 * p ’ m 2 + m 2 


0.14. 


( 6 ) 

( 7 ) 


B. The boosted r* limit 


The fact that the taus from Higgs boson decay are 
highly boosted can be exploited to obtain an even simpler 
expression for 0. We define two lightlike 4-vectors n x and 
n 2 whose spatial components, n± and fi 2 , are unit vectors 
taken along the 3-momenta of the p ± from the r decays: 

ni = (l,h p i), n 2 = (1, fip 2 ) ■ (8) 

We also define 


n\ = (1, —fipi ), n 2 = (1, —n p2 ). 


(9) 


Then, an arbitrary 4-momentum a^ can be expanded in 
terms of n x and n x as 


= 


= (a • Til) -rf- + ( a • ™l)TT + a 


-LjU 


( 10 ) 


= a + - 


"-1 




1 For a fully SU(2) X U(l) gauge-invariant form of the interac- 

tion 0 as well as an example of renormalizable realizations of 

it, see Ref. [I]. 
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where a - L by definition satisfies a 2- ■ n\ = a 2 " • n \ = 0. 
Alternatively, can be expanded in a similar manner in 
terms of n 2 and 77 , 2 . To avoid notational clutter, we do 
not make explicit whether a ± and a 2- are defined in the 
ni-rii basis or the n 2 -fi 2 basis, but it should be clear that 
the momenta of all particles originating from the r + will 
be expanded in the ni-fii basis, while those from r - will 
be in n 2 -fi 2 . 

Next, we introduce a small expansion parameter AC 1 
to quantify the degree of collinearity between the and 
the respective p . First, by construction, “collinear” 
means that p M is dominated by the + component, p + , 
while the remaining components p 2 - and p~ are small. 
Therefore, by definition, we assign the A scaling law 
p 1 - ~ Ap + to the _L component. This implies that numer¬ 
ically we have A ~ m T /mh ~ 10 -2 . Then, in the ultra- 
relativistic limit, we have 0 ~ p • p = p + p~ + p 2- • p 2 " ~ 
p + p~ + 0(A 2 p + p + ), so we see that p~ ~ A 2 p + . To sum¬ 
marize, we characterize the collinearity of p as 

(P + ,P~,P ± ) ~ (l,A 2 ,A)p+. (11) 


C. The reconstruction of neutrinos 


Since neutrinos are invisible at the LHC, we must de¬ 
vise a method to assign momenta to the neutrinos from 
the r ± decays. However, the choice of method is a deli¬ 
cate issue for the reconstruction of the 0 variable. Per¬ 
forming the A expansion in the expressions Q, we see 
that all C>(A 0 ) and ©(A 1 ) terms exactly cancel, so we are 
left with 

(fcl-p T 2) (k 2 -Pr l) - (Pt1-Pt 2) (Afi'M ~ 0( A 2 ) , 
e^upa K Pr 1 k 2 Pt2 ~ A 2 ) . 

Therefore, for the purpose of reconstructing 0, values 
must be carefully assigned to the invisible neutrino mo¬ 
menta so as not to disturb these delicate cancellations. 
In particular, the A scaling law must be respected. 

A simple method with this property is provided by the 
collinear approximation [5], where we set 

Prt = Pri = 0 > Pti =Pti = 0 , (13) 


for each i. In other words, the collinear approximation 
applies the A —>• 0 limit to the r and v momenta. Then, 
to 0( A 2 ), the expressions Q become 

cos© oc (kt~ 7 x 2 ) ( k i■ n 1 ) — (ni- 712 ) (kt-kt ), 

_L j_ ( 14 ) 

sin 0 oc €p Upc r k i p n\ k 2 p n 2 , 


again with a common proportionality factor. Finally, 
noticing that the definitions of fcf in (I5J) give fc 2 ^ = yiqi P 
in the collinear approximation, we arrive at the following 
very simple expressions that determine 0: 

cos 0 oc (qt ■ n 2 ) {q 2 ■ rii) - (m • n 2 ) (qt ■ q 2 ), 

■ lu v ±0 a l 25 -* 

sin 0 oc t pvpa q 1 q 2 F n 2 , 


with a common proportionality factor. The 4-vectors q p . 
n p and n 2 are defined in ([6]) and Q, respectively. In the 
rest of the paper, we use the expressions (15) to compute 
0 unless noted otherwise. 

One might think we should be able to do better than 
the collinear approximation by using statistical correla¬ 
tions of the neutrino momenta with other, visible, mo¬ 
menta in the process. However, the delicate cancellation 
shown in (12) implies that the correlations must get the 
neutrino momenta correctly at the level of 0( A 2 ) ~ 10 -4 
in order to outperform the collinear approximation. Al¬ 
though we could not prove that this is impossible, it cer¬ 
tainly suggests that it is extremely difficult. 


III. PARTON-LEVEL SIMULATION 

In this section, we outline our generation of the par- 
ton level samples used in the detector simulation. We 
generated tree level event samples with a combination of 
MadGraph5 (23] and MCFM as explained below for the back¬ 
ground (Z + j ) and the signal (h + j) with pTmin = 70 
GeV for the jet. The associated jet facilitates trigger¬ 
ing and provides sufficient t + t~ transverse momentum 
for the collinear approximation to hold. Other associ¬ 
ated production modes could be used, notably h + Z and 
h + W, but they have lower cross sections than h + j, 
though potentially better signal to noise that could be 
exploited in a further study. Multijet events are a poten¬ 
tial background; however, this background is found to be 
small once tau selection criteria are imposed. Therefore, 
consideration of Z +j alone is sufficient for our purposes. 

The jet-px cut was imposed at the generator level in 
anticipation of the high-px jet requirement we impose in 
our analysis to ensure that the event will pass the trigger 
as well as to avoid kinematic configurations where the 
t + t~ is produced back-to-back in the transverse plane 
and the collinear approximation for the neutrino mo¬ 
menta fails. 

The model files were generated using FeynRules pack¬ 
age Plj by modifying the SM template provided by 
FeynRules to include the CP violating coupling in the 
h-T-r vertex as well as the effective vertices for the decay 
of r —>• p + v and p —> ir + n 0 given by 

£ e ff = -p* p P^PlT - p^^dpir* - 7r*5^7r°) + c.c., 

(16) 

where P L = (1— 7s)/2 is the projection operator onto left- 
handed chirality. We implement the Higgs boson produc¬ 
tion from gg fusion using the effective field theory (EFT) 
(in which the top quark is integrated out) available from 
the FeynRules repository. In principle, such an effec¬ 
tive vertex is applicable only when the energy scale of 
the process lies well below the mass of the top quark. 
However, using the MCFM program [253, we have verified 
that the Higgs boson px distribution from the EFT in 
the region of interest (70 < px < 200 GeV) agrees at 
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the percent level with the full 1-loop QCD calculation. 
Finally, note that all the coupling constants in Ecp (16) 
are set to unity. To get the physical cross section, we use 
the cross sections from MCFM re-weighted by r -> p + v 
and p —> 7r + 7r° branching fractions. The CTEQ611 PDF 
set was used throughout. 

Parton showering and hadronization effects are in¬ 
cluded by processing the Madgraph.5 samples using the 
Pythia8 program ]26j . In principle, higher order correc¬ 
tions can be implemented by including events with higher 
jet multiplicity from tree level matrix elements and per¬ 
forming the corresponding jet merging with Pythia’s par- 
ton showering. However, we have verified that the inclu¬ 
sion of these corrections do not affect the Higgs boson px 
distribution and they differ at most by a few percent in 
the region of interest. 


IV. DETECTOR SIMULATION 

In order to incorporate detector effects into the mea¬ 
surement of the 0 angle, we use the Snowmass detector 
model [27] . implemented in Delphes 3 [28] . In Delphes, 
energy from particles deposited within the calorimeter 
is partitioned between the electromagnetic and hadronic 
sections according to a fixed ratio. This ratio is unity 
for the fraction of energy from a charged pion deposited 
in the hadronic calorimeter, which fails to account for 
energy depositions from hadronic showers that begin be¬ 
fore the hadronic calorimeter. We account for this ef¬ 
fect, including potential systematic effects due to overlap 
between the energy deposition of the neutral pion and 
the charged pion, by distributing the energy of charged 
pions in both calorimeter sections according to a proba¬ 
bility distribution that models the one given in Ref. [29] 
for 30GeV pions. 85% percent of the time the parti¬ 
cle is presumed to deposit energy via ionization only, 
and the energy is entirely deposited within the hadronic 
calorimeter. For the remaining 15% of events, the frac¬ 
tion of the particle energy deposited in the electromag¬ 
netic calorimeter is drawn from a Gaussian distribution 
of mean 0.4 and width 0.3 (selecting only physical val¬ 
ues). The rest of the energy is assigned to the hadronic 
calorimeter. The Gaussian models the probability with 
which the pion interacts early, in which case a significant 
fraction of its energy will be deposited in the ECAL. 

Tau identification in Delphes is abstracted to an ef¬ 
ficiency, which is insufficient for our purposes. There¬ 
fore, we develop an algorithm to identify taus explicitly 
that assigns the charged pion and the neutral pion to a 
track and a calorimeter tower, respectively. We select 
jets reconstructed with the anti-fcx algorithm [30] with 
a distance parameter of 0.5, which are reconstructed in 
Delphes using FastJet m ■ At least three jets are re¬ 
quired in each event with px > 30GeV within |p| < 2.5 
(the tracking limit of the detector), of which the lead 
jet must have px > 100 GeV (well above the parton level 
requirement of 70 GeV). All jets with px >30 GeV are ex¬ 


amined as potential r-candidates. The lead track within 
A R < 0.2 of the jet axis is selected, and the sum of 
the transverse momenta of the remaining tracks, within 
A R < 0.4, is required to be less than 7 GeV in order to 
exclude jets that are unlikely to be taus. These selection 
requirements emulate the more sophisticated selection of 
taus in the LHC experiments, and are treated as equiva¬ 
lent to the Delphes efficiency of 65%. 

We make further selections on the individual elements 
within the tau candidate in order to ensure that the decay 
products have been measured with sufficient precision for 
the calculation of 0. The lead track is required to have 
px > 5 GeV, or this candidate is rejected. We then pro¬ 
ceed to examine calorimeter towers within A R < 0.2 of 
the jet axis. We select the tower with the largest trans¬ 
verse energy in the electromagnetic portion of the tower. 
The tower in question must also have an electromagnetic 
fraction (EMF) of 0.9, that is, at least 90% of its total 
energy within the ECAL. If this tower has px > 5 GeV, is 
within A R < 0.2 of the selected track, and the sum of the 
selected tower transverse energy and the track momen¬ 
tum is greater than 20 GeV, then this jet becomes a r- 
candidate suitable for 0 reconstruction, and the charged 
pion is assumed to be the selected track with the neutral 
pion the selected tower. 

Finally, the electromagnetic portion of the actual AT¬ 
LAS and CMS detectors have finer granularity than a 
single tower, which is not properly accounted for by the 
Snowmass detector. To approximate a position resolu¬ 
tion of one detector element of the LHC detectors, we 
introduce the resolution for the position in p-f/ of the 
neutral pion as a Gaussian smearing with standard de¬ 
viation 0.025/\/l2. A resolution no better than a tower 
(as in the default Snowmass detector) was found to de¬ 
grade significantly the sensitivity to the modulation in 
the 0 distribution. Finally, events are required to have 
at least two r-candidates with opposite charges. Should 
there be more than two r-candidates, the highest px pair 
with opposite charge is selected. 


V. RESULTS 

We now come to the critical question of whether 
the signal modulation survives detector effects. In or¬ 
der to answer that question, consider the analysis cut- 
flow shown in Table [I] The signal acceptance - the ra¬ 
tio of the second number to the first in column two 
of Table |T] - is 0.430, while that of the background is 
0.254. The products of the efficiencies of all the re¬ 
maining cuts including the 65% r reconstruction effi¬ 
ciency are 0.65 2 x 143979/1078279 = 0.0564 and 0.65 2 x 
134628/1069747 = 0.0531 for the signal and background, 
respectively. Therefore, overall, the signal and back¬ 
ground efficiencies are 0.430 x 0.0564 = 2.42 x 10 -2 and 
0.254 x 0.0531 = 1.35 x 10 —2 , respectively. 

To see the effects of finite detector resolution on our 
ability to reconstruct the signal 0 distribution ([3]), we 
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Selection criterion 

Z+j-> TT+j 

ft+j-T rr+j 

None 

9618732 

4698651 

At least three jets 

Pt > 30 GeV, lead jet 
Pt > 100 GeV 

2445860 

2019316 

Two r-candidates 
only AR and 
isolation requirements 

1069747 

1078279 

r-track pt > 5 GeV 

903429 

941368 

Tower Et > 5 GeV 

804566 

856299 

Tower EMF >0.9 

156140 

160312 

Sum pt > 20 GeV 

143508 

149591 

Opposite Charge 

134628 

143979 


TABLE I: Cut-flows for Z —> tt and h — > tt. The generator 
and reconstruction level cuts on the jet are pt > 70 GeV, and 
100 GeV, respectively. See text for more details. 



FIG. 1: Comparison of parton level 0 distribution (Monte 
Carlo Truth) with the distribution after the Delphes simula¬ 
tion and event reconstruction for h — ► tt. 


parameters 

MC Truth 

Delphes Sim. 

c 

0.0397 ±0.0001 

0.0397 ±0.0002 

a 

0.0089 ±0.0001 

0.0085 ± 0.0002 

7 

1.03 ±0.01 

1.03 ±0.02 

d 

-0.003 ±0.003 

0.01 ±0.01 

a = a/c 

0.225 ±0.002 

0.214 ±0.004 


TABLE II: Fits to the signal 0 distribution before and after 
the inclusion of detector effects using the functional form c — 
a cos( 7 © — 2<5). 


present in Fig. [T] the reconstructed 0 distribution from 
the signal samples generated with A = 0 before and after 
the inclusion of detector effects, labelled as Monte Carlo 
Truth and Delphes Simulation , respectively. We perform 
fits to these distributions using the functional form c — 
a cos( 7 @ — 2(5) with free parameters a, c, 7 , and S in 
order to extract the modulation amplitude a = a/c. The 
fit results are presented in Table [TT] and show that 7 ~ 1 
and S ~ 0 as expected from the analytical result © and 
the assumption A = 0, even after detector effects. Most 
crucially, the dilution of the modulation amplitude a due 
to the detector effects is only about ~ 4%, going down 
from 0.225 to 0.214. 

In order to check that our r reconstruction algorithm 
does not introduce artificial modulations in the back¬ 
ground 0 distribution, we present in Fig. [2] the Monte 
Carlo Truth and Delphes Simulation 0 distributions us¬ 
ing the background samples. Clearly, the background is 
consistent with a flat distribution after the detector ef¬ 
fects and we may conclude that the latter do not bias the 
shape of the 0 distribution. 

Now, we address the crucial question: how well can we 
distinguish the 0 distribution with A = 0 from that with 
A 7 ^ 0? Given N events, each yielding a measurement of 
0 ,:, we construct a test statistic t designed to distinguish 
between the SM hypothesis A = 0 and the alternatives 



FIG. 2: Comparison of parton level 0 distribution (Monte 
Carlo Truth) with the distribution after the Delphes simula¬ 
tion and event reconstruction for Z — >• tt. 


A 7 ^ 0. A standard choice for t is 


t = 


2 In 


L(A) 

m 


where 


N 

L(A) = n/(0<|A) 

i=1 


(17) 

(18) 


is the likelihood function, while /(0j|A) is the 0 distri¬ 
bution for a given value of A. 

In order to enhance the signal to background ratio 
( S/B ), we use a boosted decision tree (BDT), trained 
using the invariant mass of the reconstructed r ± system, 
the kinematic variables (px, (j>, v) °f the final state visible 
particles, and the missing transverse energy. Therefore, 
each event can be characterized by the measurement of 
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©i and the value of the BDT discriminant d,, which yields 
the modified likelihood function, 

N 

L{A) = II /(©i, d»|A), (19) 

i— 1 

where the probability density function of the data 
{(©, d)} is modelled as follows, 

/(©, d|A) 

V ' (20) 

oc B Pbdt(^) + S Pbdt(^) [l - a cos(0 - 2A)] , 

where Pbdt(^) an d Pbdt(^) are > respectively, the signal 
and background BDT distributions for the discriminant 
parameter d, while B and S are, respectively, the total 
expected background and signal counts given by 



B = e B a B C, S = es<JsB. 


( 21 ) 


Here, as and a B are the associated cross section times 
branching fractions, the e’s are the signal and background 
efficiencies discussed above and C is the integrated lumi¬ 
nosity. 

In order to quantify how well one might be expected 
to distinguish a A/0 hypothesis from the SM hypoth¬ 
esis A = 0, we generate three samples of 0 and d val¬ 
ues with A = 0, 7 t/ 4, and n/2. For each sample, we 
calculate the values of t as defined in Eq. (17), which 
yield the t distributions as illustrated in Fig. |3| for the 
A = 0 and 7 t/ 4 samples. To quantify the discrimination 
between the two hypothesis, the p-value of the SM t dis¬ 
tribution is computed using the median of the non-SM 
t distribution. As proxies for systematic uncertainties in 
the S and a parameters, the p-values are varied in the 
10% neighborhood of their nominal values. In Fig. [4] we 
present the reach as a function of the integrated lumi¬ 
nosity for A = 7 t/ 4 and tt/2. Because we have chosen 
the “observed” value of t to be the median of the non- 
SM hypothesis, the reach in this context is interpreted 
as follows: if the non-SM hypothesis is true, then there 
is a 50% chance to reject the SM hypothesis A = 0 at 
the Z-sigma level, where Z = d*” 1 ^) = v / 2erf^ 1 (2 p— 1) 
with p being the p- value shown in Fig. [3] 

It is clear from Fig.[4]that the optimistic conclusion ob¬ 
tained for an ideal detector does not survive the smearing 
effects of a realistic detector and our reconstruction algo¬ 
rithms. At best, at ~ 1500 fb — 1 , there is a 50:50 chance 
of rejecting the A = 0 hypothesis at 95% confidence level. 
It is important to note where the poor performance stems 
from. As mentioned earlier, the reduction of the modu¬ 
lation amplitude a of the signal due to the angular res¬ 
olution is negligible. Rather, we find that the discrim¬ 
ination power of the BDT dominantly comes from the 
invariant mass of the reconstructed t + -t~ system, and 
consequently the S/B ratio is sensitive to the mass reso¬ 
lution, which in turn is sensitive to the resolution of the 
missing transverse energy (MET). Therefore, our anal¬ 
ysis shows that the current MET resolution of the LHC 


FIG. 3: Distributions of the t statistic for the SM hypothesis 
(A = 0) and a non-SM hypothesis (A = 7r/4). The “ob¬ 
served” value of the statistic is taken to be the median of the 
t distribution of the non-SM hypothesis. 



FIG. 4: The reach, quantified on the Z-sigma scale of the 
normal distribution, versus the integrated luminosity. The 
uppermost plot shows the effect of reducing the background 
by a factor of two, while keeping every else fixed. 


detectors is too low to achieve a sufficient reduction of 
the total background count B. 

To illustrate the effects of improved missing transverse 
energy resolution, we show in the uppermost plot in Fig. [4] 
the effect of increasing S/B by a factor of two by decreas¬ 
ing B by a factor two, while keeping everything else the 
same. Increasing the S/B by increasing S' by a factor 
two (while keeping a unchanged) has a similar but more 
pronounced effect. 
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VI. CONCLUSIONS 

We have examined a recent proposal [1] for measur¬ 
ing a potentially non-zero CP phase A in Higgs to r + r” 
decays using an angular variable 0 that can be recon¬ 
structed from the decay chains r ± —>• p ± b l 11 —> n ± TT 0 i i . A 
non-zero phase of the h-T-r Yukawa coupling would be 
unambiguous evidence of physics beyond the SM. Con¬ 
trary to previous expectations in. with ~ 1000 fb 1 we 
find that the standard model CP phase, A = 0°, can be 
distinguished from the A = 90° hypothesis at no more 
than 95% confidence level. 

We have addressed two obvious effects that are ex¬ 
pected to degrade our ability to measure A from the 0 
distribution. First, we examined the issue of neutrino 
momentum reconstruction, for which the original pro¬ 
posal [T] employed the collinear approximation. We have 
presented a theoretical argument that strongly suggests 
that this cannot be improved further. 

Second, we have closely studied how finite resolutions 
of a realistic detector affects our ability to distinguish be¬ 
tween the SM signal and its deviations induced by a non¬ 
zero CP violating phase A, using the Snowmass detector 
model implemented in the Delphes detector simulation 
program. Contrary to the naive expectation that a finite 
angular resolution should hamper the reconstruction of 
0 as it requires resolving the charged and neutral pions 


from the r decays, we have found that the finite angular 
resolution of the LHC detectors reduces the modulation 
amplitude of the signal 0 distribution by only ~ 4%, 
compared to that of an ideal detector. Instead, we have 
found that the experimental ability to distinguish the 
hypotheses A = 0 from A / 0 is strongly limited by 
the missing transverse energy resolution. The finite en¬ 
ergy resolution significantly degrades the use of the rr 
invariant mass distribution as the principal discriminant 
between the h —» rr signal and the Z —»• rr background. 
In contrast, the reconstruction of 0 itself is nearly unaf¬ 
fected by the finite energy resolution, as the 0 variable 
is basically determined by the angular topology of the 
Higgs boson decay products. Fortunately, the resolution 
of missing transverse energy is an area in which improve¬ 
ments can be hoped for, e.g., by using multivariate meth¬ 
ods |32]. 
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